arXiv:1509.03526v2 [hep-ex] 14 Dec 2015 


EUROPEAN ORGANIZATION EOR NUCLEAR RESEARCH 



CERN-PH-EP-2015-245 


September 11, 2015 
rev. Dec. 14, 2015 


Longitudinal double spin asymmetries in single hadron 
quasi-real photoproduction at high pr 


Abstract 

We measured the longitudinal double spin asymmetries An for single hadron muoproduction off 
protons and deuterons at photon virtuality <\ (GeV/c)^ for transverse hadron momenta pj in the 
range 1 GeV/c to 4 GeV/c. They were determined using COMPASS data taken with a polarised muon 
beam of 160 GeV/c or 200 GeV/c impinging on polarised ®LiD or NH 3 targets. The experimental 
asymmetries are compared to next-to-leading order pQCD calculations, and are sensitive to the gluon 
polarisation AG inside the nucleon in the range of the nucleon momentum fraction carried by gluons 
0.05 <Xg < 0.2. 
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1 Introduction 

The spin structure of the nucleon is one of the major unresolved issues in hadronic physics. While the 
quark spin contribution to the nucleon spin, denoted as AT, has been measured to be about 30% lUl, the 
gluon spin contribution is still insufficiently constrained after more than two decades of intense study. In 
the framework of perturbative Quantum Chromodynamics (pQCD), inclusive Deep Inelastic Scattering 
(DIS) is sensitive to gluon contributions only through higher-order corrections to the cross section. The 
spin-averaged gluon density g(xg), where Xg denotes the nucleon momentum fraction carried by gluons, 
is well constrained by DIS experiments with unpolarised beam and target because of their high statistics 
and large kinematic coverage. The fewer data from DIS experiments with polarised beam and target, 
however, can not sufficiently constrain the gluon helicity distribution Ag(xg). This affects directly our 
knowledge of the contribution of the gluon spin to the spin of the nucleon, known as AG = / Ag(xg)dxg, 
and to a lesser extent that of the quarks |2S- In order to better constrain Ag(xg), one has to resort to 
processes where contributions from gluons appear at leading order, such as hadron production at high 
transverse momenta or production of open charm in polarised lepton-nucleon 131IH |5l |6l |71 or hadron- 
hadron interactions ll^l9l lT0llTn . 


The COMPASS collaboration has already investigated asymmetries of hadrons at high transverse mo¬ 
menta pt, in both the DIS and the quasi-real photoproduction regimes ||4l|6l|T2l. Here, transverse means 
transverse with respect to the direction of the virtual photon y* that is exchanged in the scattering process. 
Using a Lund Monte Carlo simulation, these measurements were interpreted on the hadron level, thereby 
simultaneously extracting the gluon helicity on the parton level. Such an analysis is restricted to leading 
order (LO) in the strong coupling constant a^, as presently there exists no next-to-leading order (NEC) 
Monte Carlo simulation for leptoproduction. Due to the limitation of neglecting gluon contributions at 
NEC, such results can not be used in recent global hts at NEC of polarised Parton Distribution Eunctions 
(PDEs) (H. 

In this Letter, we present a new analysis of COMPASS data for single-inclusive hadron quasi-real pho¬ 
toproduction at high Pt which differs from our previous analysis in that all measured hadrons within a 
given Pt bin are included in the analysis, and not only the hadron(s) with highest pT- Moreover, the inter¬ 
pretation of the results is based on a collinear pQCD framework that was developed up to NLO llT4l . the 
basic concept being the application of the factorisation theorem to calculate the cross section of single- 
inclusive hadron production. The authors of Ref. ifHIl discuss the sensitivity of COMPASS data to Ag{xg) 
in terms of contributions from “direct-photon”, y^g —)• qq (Photon Gluon Eusion), and from “resolved- 
photon” subprocesses, qg and gg, where the photon acts as a source of partons. Similarly, they consider 
direct y*q —)■ qg (QCD Compton) as well as resolved qq and gq subprocesses for the background. These 
contributions to the cross section are represented schematically in Eig. [T] In the framework of collinear 
fragmentation, photo-absorption on quarks, 7 *^ —is not contributing to high pT hadron production. 


'Note that also inclusive quasi-real photoproduction of hadron pairs can be considered HU. 
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Fig. 1 : Contributions to the single-inclusive cross section for quasi-real photoproduction of a hadron h into direct 
(left) and resolved (right) subprocesses according to Ref. 03. The internal lines represent the photon a — y* (left) 
and partons {a,b,c} = {q,q,g}- The central blob describes the hard scattering cross section a. The peripheral 
blobs describe the non-perturbative objects; parton distributions of the nucleon, and of the photon, /J*, and 
the fragmentation functions of the produced hadron, D*. 


In order to gain confidence in the applicability of this pQCD framework to single-hadron production 
with longitudinally polarised beam and target, an important step is to compare predictions of this model 
to measurements with beam and target unpolarised, for which the PDFs are well known. While good 
agreement was found by RHIC experiments on the production of high-pr hadrons in pp collisions 
at ~ 200 GeV ifT^ [TTll . complications arise when hard scattering subprocesses are probed in the 
“threshold” regime, in which large logarithmic contributions from soft and collinear gluons play a sig¬ 
nificant role ifTSl . Such contributions become dominant at the COMPASS centre-of-mass energy of 
~ 18 GeV. When taken into account by a technique known as “threshold resummation” at next- 
to-leading logarithm (NLL) lIT^ . the calculations reproduce the COMPASS cross section measure¬ 
ments lfT9ll within theoretical uncertainties. 

In this Letter, we analyse the quasi-real photoproduction data collected by COMPASS from 2002 to 
2011 on longitudinally polarised deuteron and proton targets. In Sec.|^ we give a brief description of 
the experimental setup, and details on the data selection can be found in Sec. The procedure for the 
asymmetry determination is described in Sec. In Sec. we present the corresponding double spin 
asymmetries for single-inclusive hadron production as a function of their transverse momenta pj. These 
asymmetries are compared to calculations that were performed using the code of Ref. Ifldl . which does 
not include the resummation of threshold logarithms. 


2 Experimental Setup 

The measurements were performed with the COMPASS setup using positive muons from the M2 beam 
line of the CERN SPS. A detailed description of the experimental setup can be found in Ref. Il20ll . 
with updates valid since 2006 described in Ref. |[2ll . The muon beam had a nominal momentum of 
160 GeV/c, except for 2011 where the momentum was 200 GeV/c. On average, its momentum spread 
was 5 % and its polarisation was ss 0.8. Momentum and trajectory of incident muons were measured 
by a set of scintillator hodoscopes, scintillating fibre and silicon microstrip detectors. The beam was 
scattered off a solid state deuterated lithium (®LiD) target from 2002 to 2006 and off an ammonia (NH 3 ) 
target in 2007 and 2011, providing longitudinally polarised deuterons and protons, respectively. The 
target material was placed inside a large aperture superconducting solenoid, and by dynamic nuclear 
polarisation it was polarised to a value of P 0.5 for ^LiD and P 0.85 for NH 3 . Until 2004, the 
target material was contained in two contiguous 60 cm long cells that were oppositely polarised. From 
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2006 onwards, three contiguous target cells of length 30cm, 60cm and 30cm were used to minimise 
systematic effects, with the polarisation in the outer cells being opposite to that in the central one. The 
direction of the target polarisation was regularly flipped by reversing the solenoid field to compensate 
for acceptance differences between the different target cells. At least once per year, the direction of the 
polarisation was reversed relative to that of the solenoid field. The dilution factor /, which accounts for 
the presence of unpolarisable material, amounts to typically 0.4 for the deuterated lithium target and to 
0.16 for the ammonia one. It is calculated as the ratio of the cross section on polarisable deuteron or 
proton to that on all target nuclei, corrected for unpolarised x and y dependent electromagnetic radiative 
effects ll22]l . where .r is the Bjorken scaling variable and y the relative muon energy transfer. No further 
radiative effects are taken into account. The momenta and angles of scattered muons and produced 
hadrons were measured in the two-stage open forward spectrometer, where each stage includes a dipole 
magnet with upstream and downstream tracking detectors. 


3 Data Selection 

In order to be selected, an event must have an interaction vertex that contains both incoming and scattered 
muons and at least one hadron candidate track. The measured beam momentum is required to be in a 
±20 GeV interval around the nominal value (±15 GeV in 2011). In order to equalise the flux through 
each target cell, the extrapolated beam track is required to pass all target cells. Cuts on the position 
of the vertex allow the selection of the target cell, in which the scattering occurred. Only events with 
photon virtuality < 1 (GeV/c)^ are accepted. This kinematic region is referred to in the following 
as quasi-real photoproduction region. In addition, y is required to be within 0.1 and 0.9, where the 
lower limit removes events that are difficult to reconstruct and the upper limit removes the region where 
electromagnetic radiative effects are large. These kinematic cuts result in a range of 10^^ < x < 0.02 and 
a minimum mass squared of the hadronic final state, W^, of 25 (GeV/c^)^. The hadron candidate track 
must have pr > 0.7 GeV/c. The fraction z of the virtual photon energy carried by the hadron is required 
to be in the range 0.2 < z < 0 . 8 , where the lower limit is imposed to suppress the contribution from target 
remnant hadronisation and the upper limit to reject badly reconstructed hadrons. The angle between the 
direction of the hadron and that of the virtual photon is restricted to be in the range lOmrad < 6 < 
120mrad, which corresponds to 2.4 > rj > —0.1, where T] is the pseudo-rapidity in the 7*A centre-of- 
mass system. After all selections, the final sample consists of 140 million events for the deuteron target 
and 105 million for the proton target. 


4 Asymmetry Calculation 

The double-spin asymmetry of the cross sections for single hadron quasi-real photoproduction is de¬ 
fined as All = ± = Aa/a, where the symbols ^ and denote anti-parallel 

and parallel spin directions, respectively, of the incident muon and the target deuteron or proton. This 
asymmetry is evaluated using the same method as in our previous analyses JQ. The number of hadrons 
produced in a target cell is related to All and to the spin independent cross section a = ± G^: 

Ni = a;0;n;a(l + fiPbPtALh), where i = u\,di,U 2 ,d 2 - A target cell (u or d) with a given direction of the 
target polarisation (1 or 2 ) has the acceptance a;, the incoming muon flux 0 ,- and the number of target 
nucleons n,. For the two-cell target, u and d denotes upstream and downstream cell, respectively, while 
for the three-cell target, u denotes the sum of the outer cells and d the central cell. The asymmetry All 
is extracted from the second order equation that is obtained from the quantity ■Nd 2 )/{Ndi -Auj). In 
this relation, fluxes and acceptances cancel, provided that the ratio of acceptances of the two sets of cells 
is equal for the two orientations of the solenoid field. 

In order to minimise statistical uncertainties, all quantities entering the asymmetry are calculated for 
each hadron using a weight factor w, = fiPt (231. The muon beam polarisation Ph is obtained from a 
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parametrisation as a function of the beam momentum. The target polarisation Pf. is not included in the 
weight Wi as it changes with time and could generate false asymmetries. In order to reduce systematic 
uncertainties, data are grouped into periods that are close in time and hence have the same detector condi¬ 
tions, and the weighted average over all periods is taken. The asymmetries determined for a given target 
from data taken in different years were found to be consistent and hence combined. The asymmetries 
are obtained for both positive and negative unidentified hadrons in bins of pj in the range 0.7 GeV/c 
to 4 GeV/c and in bins of T] in the range —0.1 to 2.4 in order to facilitate a detailed comparison to the¬ 
ory (see Sec. [^. The data for pj- < 1.0 GeV/c are only used to investigate systematic uncertainties as 
the pQCD framework is commonly applied only for hard scales ^ p\'> 1.0(GeV/c)^, and they are 
shown greyed out in all the figures where they appear in. 

The systematic uncertainties on An are calculated as the square root of the sum of squares of multiplica¬ 
tive and additive contributions. The uncertainties on the dilution factor (?a 5%), the beam (« 5%) and the 
target (ss 5%) polarisations contribute to a total of ss 8% of multiplicative uncertainties, i.e. those being 
proportional to the asymmetry itself. Additive contributions originate from fluctuations of the detector 
performance, which may lead to false asymmetries. Their possible occurrence is investigated by dividing 
the data sample into different subsets. Asymmetries calculated with hadrons detected in left and right 
(top and bottom) parts of the spectrometers are found to be compatible within statistical uncertainties, 
as well as those for the two relative orientations of the solenoid field and the target spin vectors. No 
systematic uncertainty is thus attributed to these effects. Possible false asymmetries between data sets 
having the same polarisation states are also found to be compatible with zero. For each pj bin, the sta¬ 
tistical distribution of the asymmetries calculated by time periods closely follows a normal distribution. 
The observed deviations from a Gaussian allow us to quantify the level of overall additive systematic 
uncertainties as a fraction of the statistical ones, which on average amounts to about one half. These 
additive systematic uncertainties largely dominate over the multiplicative ones. 

5 Results and interpretation 

The final asymmetries are calculated using all data accumulated with the deuteron target in the years 
2002 to 2006 and with the proton target in the years 2007 and 2011. Their pr-dependence in three 
rapidity bins spanning the full interval —0.1 < rj < 2.4 ([—0.1,0.45], [0.45,0.9], and [0.9,2.4]) is shown 
in Fig. and Fig.[^for each target type and hadron charge. 

We compare our asymmetries with theoretical calculations at NLO without threshold resummation based 
on the framework described in Ref. |[T4l and summarised in the following. Using the code of Ref. |[T4ll . 
the asymmetries are computed in bins of pr and rj as the ratio of polarised to unpolarised hadron cross 
sections, where a cross section is a convolution of the “muon-parton distribution function” fa, the nu¬ 
cleon PDFs f^, the perturbative partonic cross sections da+b->-c+x> and the fragmentation functions (FF) 
£)/!. 


All[pt,^) 


dAa'’ 

da* 


{PT,r]) 


'La,b,c ^fb dAdg+h^c+X Dc 

'La.bxfa ® fb ®<^^a+b^c+X®D’^c 


( 1 ) 


Here and below, spin-dependent quantities are denoted by the symbol A and will be referred to as po¬ 
larised ones in the rest of the Letter (spin-independent ones as unpolarised). The processes involved in 
Eq. [T]can be classified into “direct” ones that are initiated by a quasi-real photon and “resolved” ones 
that are initiated by its fluctuation into partons. This classification is denoted by the subscript a (see 
Fig. j^. For direct processes, subscript a refers to y*, and (A)/^» is the probability for a muon to emit a 
quasi-real photon. For resolved processes, subscript a refers to q, q or g, and {A) fa is the convolution of 
this probability with a non-perturbative parton distribution of the photon, {A)f^ . The polarised version 
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Fig. 2: The asymmetry An as a function of pj for charged hadron quasi-real photoproduction on the proton for 
three rapidity bins. The bands at the bottom indicate the systematic uncertainties, which are dominated by time 
dependent fluctuations (see Sec.[^) Top: positive hadron production; bottom: negative hadron production. 














































of the latter is not known experimentally and hence taken to range between the positive and negative 
magnitude of the unpolarised one. This induces a small uncertainty in the theoretical calculations. 


The values of the asymmetries are computed here using the following input distributions: the unpo¬ 
larised parton distributions of the photon from GRS ll2^ . the unpolarised nucleon PDFs from 
CTE065 l[^ . the three polarised PDF sets from GRSV ll26l as in Ref. lIT^ (the “standard” set and the 
two sets for “maximum” [Ag(x) = g(x)] and “minimum” [Ag(x) = —g(x)] gluon distribution functions 
at input scale), as well as the most recent polarised PDF set DSSV14 from Ref. llT3]l . For the polarised 
PDF sets used, the integration over the range 0.05 <Xg < 0.2, which is characteristic for the kinematic 
coverage of COMPASS in the gluon momentum fraction, yields the following “truncated” values of AG 
at a scale of 3 (GeV/c)^ : AGorsv^,. ~ —0.6, ^Gdssvu ~ 0.1, AGgrsv„j ~ 0.2, AGgrsv^oj, ~ 0.7. The 
other inputs that we changed with respect to Ref. ifTdll are the fragmentation functions for which we 
use the most recent parton-to-pion fragmentation set of Ref. |[27l . which best fits the recent COMPASS 
pion multiplicities ||281 . We checked, as it was done in Ref. |[T4ll . that asymmetries for hadron and pion 
production are almost indistinguishable, so that it is safe to compare our experimental data to theoretical 
asymmetries computed with a parton-to-pion FF set. The fractions of the total unpolarised (respectively 
polarised) cross section for the various individual subprocesses are shown in Fig. as a function of pr- 
Although the unpolarised cross sections for processes involving gluons from the nucleon are not the dom¬ 
inant ones, the polarised cross section for the yg subprocess is large in magnitude for hadron production 
at high pt, which makes the study of such asymmetries relevant in the COMPASS kinematic region. 
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Fig. 4: Contributions of the six subprocesses a + b ^ c + X described in Sec. [^to the full NLO unpolarised 
(left) and polarised (right) photoproduction cross sections for a deuteron target. The polarised cross sections are 
computed using the polarised PDF set of Ref. M- 


The computations of All{pt) are performed at COMPASS kinematics using the same cuts as in the 
present data analysis, i.e., pr > ^ GeV/c, < 1 (GeV/c)^, 0.1 < y < 0.9, 0.2 < z < 0.8. For consistency, 
we verified that when using the same inputs as in Ref. ITHI . we reproduce the asymmetry calculated there. 
The computations are done separately for the production of positive and negative hadrons, and for three 
distinct bins in 77 : [—0.1,0.45], [0.45,0.9], and [0.9,2.4]. The results of these computations are compared 
to the experimental asymmetries in Fig.J^and Fig.[^ 

The data are seen to be consistent with the NLO calculations of Ref. |[T4l using the most recent polarised 
PDF mi and FF lETl sets, except for positive hadron production from the proton in the rapidity range 
—0.1 < rj < 0.9. Our data are also compared in Fig. [^and Fig. to calculations using earlier GRSV 
polarised PDFs ll26l to give an impression of their sensitivity to AG, which seems enhanced at higher 
values of T]. A possible reason for the discrepancy seen at low rj in positive hadron production from 
the proton may be that the model calculations had to be done without threshold resummation at NLL, as 
the formalism for the polarised case is not yet fully available. However, contrary to the unpolarised case 
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Fig. 5: COMPASS asymmetries An for a proton target as a function of pj and in three rapidity bins, compared 
to NLO calculations based on Ref IH for different choices for the polarised PDFs (see text). Only statistical 
uncertainties are shown. Error bands on the theory curves represent the uncertainties due to the polarised parton 
distribution of the photon. Top; positive hadron production; bottom; negative hadron production. 
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Fig. 6: Same as Fig.|^ but for a deuteron target. 
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where data can only be described by pQCD if threshold resummation is included |[T8]| ^ spin asymmetries 
are expected to be less affected 1291 . A first estimation of the impact of threshold resummation on our 
double spin asymmetries including so far only the direct processes ll^ indicates a substantial dilution of 
the asymmetries, which may explain part of the discrepancy between experiment and theory for positive 
hadron production on the proton at low values of rj. 

6 Summary 

In summary, we have presented in this Letter a new analysis of COMPASS data on polarised single- 
inclusive hadron quasi-real photoproduction, which in principle is well suited for an extraction of the 
gluon polarisation AG in the framework of collinear pQCD. Results for the longitudinal spin asymmetry 
AllIpt) on polarised protons and deuterons are given separately for positively and negatively charged 
hadrons, and in three rapidity bins. They are compared to theoretical calculations at NLO without thresh¬ 
old resummation and overall agreement is found with the calculations based on earlier GRSV^f^ and 
recent DSSV14 polarised PDF sets, and using the most recent FF set. Nevertheless, calculations includ¬ 
ing full threshold resummation at NLL are needed before a meaningful result on AG can be extracted 
quantitatively from our data. 
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